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We propose a simplified mathematical construction of the quantum Fourier transform which is 
suited for systems described by Ising-type Hamiltonians. By contrast to the standard Cooley-Tuckey 
scheme, which prescribes sequences of CPHASE gates, our implementation is based on one-qubit 
gates and a free evolution process. We also show how to obtain a quadratic speed-up by applying the 
conditional interactions simultaneously. Thus rather than 0{N^) our implementation time scales as 
0{N). Finally, we show a realization of our method with homogeneous microwave driven ion traps 
in a magnetic field with gradient. 
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I. INTRODUCTION 


The quantum discrete Fourier transform (QFT) is the 
fundamental ingredient of many principal quantum al¬ 
gorithms, known to provide substantial speed-up over 
their classical counterparts [1-3], including Shor’s or¬ 
der finding, factorization and the discrete logarithm [4], 
Deutsch’s algorithm [5] and phase estimation [6]. Being 
at the heart of quantum computation QFT has become 
a textbook example of the potential superiority of quan¬ 
tum over classical computers and thus it has inspired 
extensive research aiming at the implementation of this 
transformation in different physical systems. 

Experimental demonstrations include the composition 
of two- and three-qubit QFT in nuclear magnetic reso¬ 
nance (NMR) [7-10], neutral molecules [11], supercon¬ 
ducting qubits [12], and semiclassical realizations with 
trapped ions [13] and photonic qubits [14]. These real¬ 
izations use predominantly the Cooley-Tuckey algorithm 
[15], which gives an efficient implementation of QFT with 
a polynomial number 0{N‘^) of conditional gates, with N 
being the number of qubits. Despite the mathematical ef¬ 
ficiency of this algorithm, however, the attempts to scale 
the number of qubits have quickly come across a serious 
obstacle: the number of gates and the ensuing physical 
operations with limited precision still grow rapidly with 
the size of the problem. They reach prohibitively large 
values before approaching even moderate computational 
scales. Therefore one must seek for more realistic im¬ 
plementation models which require less resources. One 
way to go is by decomposing QFT directly into more 
favourable gates, which can be realized naturally in the 
particular physical system. 

In this work we present a mathematical construction 
of QFT, which is a modification of the Cooley-Tuckey 
scheme. Our circuit is suited for a broad range of systems 
described by Ising-type Hamiltonians of the form 


Hi 


h 

2 




k<.l—l 


where represent the spins and Jki represent coupling 
strengths. Examples for such systems are nuclear mag¬ 
netic resonance [8], ion traps [16], spin chains. Our pro¬ 
posal has two principal advantages. Eirst, the sequences 
of conditional gates, which are the major stumbling block 
to QFT, are now carried out by simple time evolutions 
driven by the Hamiltonian (1) over certain time intervals. 
Thus in our implementation QFT is realized with a se¬ 
ries of one-qubit gates applied at particular times. This 
implies that, second, the number of physical interactions 
with the system is reduced. We also show how to obtain 
quadratically faster QFT by applying the conditional in¬ 
teractions simultaneously. 

We include an explicit discussion of the physical real¬ 
ization of the proposed scheme with trapped ions. 


II. QUANTUM FOURIER TRANSFORM: 
MATHEMATICAL CONSTRUCTION 


The iV-dimensional QFT is defined through its action 
on the computational basis states JO), Jl), [2), ..., |iV—1) 
by 


1 

= ( 2 ) 

It is a unitary transformation, which maps each state jn) 
into an equal superposition of all states each imprinted 
with a particular phase. 

A circuit that implements QFT is shown in Fig. I 
(top). It is based on the Cooley-Tuckey algorithm and 
can be found in Refs. [I] and [17]. Two gates from the 
elementary gate set of the standard computational model 
are used - the single-qubit Hadamard gate 
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acting on qubit k and the two-qubit conditional phase 
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gate (CPHASE gate) 
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acting on qubits k and /, where 
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(5) 


This circuit offers an exponentially faster implementation 
(with 0{N‘^) gates) compared to the best known classi¬ 
cal algorithm for fast Fourier transform (with 0{N2^) 
gates). Despite this formal superiority, however, this im¬ 
plementation has proved difficult to scale due to the se¬ 
quences of the CPHASE gates. Hence other mathemat¬ 
ical constructions of QFT are desirable, which use more 
favourable gates for the particular physical system rather 
than the CPHASE gates. 

Based on the Cooley-Tuckey scheme, in the following 
we give a mathematical construction of QFT, which is 
adapted for a class of systems described by the Hamilto¬ 
nian (1). This Hamiltonian realizes long-range spin-spin 
coupling and, since it is time-independent, the evolution 
matrix \s U = exp {—iHit/H), or 


U = exp 





( 6 ) 


Below we will also use the notation U (t) = 


exp 


^^z2k<l=l '^klO'z Oz 


when we want to highlight 


the dependence of U on time. Note that U yields 
control-phase shifts between all pairs of spins, U = 

nf<.= iUki{(t)ki), where Ukii'pki) = exp . 

The corresponding phase is (j)ki = Jkit/2. Thus Uki{4>ki), 
which will substitute the CPHASE gate, is merely de¬ 
livered by free evolution after time t = 2(j)ki/Jki- Thus, 
the required free evolution time might be different for 
each (pki; however, if the coupling Jki can be tailored 
such, that all t are identical, all phases (pki can be ac¬ 
quired in a single step. In NMR, there exist proposals 
for the implementation of the two-qubit CNOT gate and 
multi-qubit gates, based on always-on interactions (by 
refocussing other interactions, and thus “wasting them”) 
[18]. 

Our gate set comprises a rotation Rk (</>), a phase gate 
Tk[4>)^ and the entangling gate Uki{4>)- The indices in¬ 
dicate qubit numbers. We will also use implementa¬ 
tions, where the phase gate is incorporated in the ro¬ 
tation to form R{0, (p) according to the rule R{0^ p) = 
T{(j)/2)R{9)T{—(f>/2). Such phased rotation can be ob¬ 
tained with a simple phase shift (p in the driving field 
implementing R{0). We have 


Rk{e,(l)) = 

(7a) 

Tk{4>) = 

(7b) 

II 

(7c) 



FIG. 1: (top) Decomposition of the quantum Fourier trans¬ 
form using standard gates from the circuit model [1, 4]. The 
example shown is for N = 4 qubits. Here H is the Hadamard 
gate and C{4>k) are control phase gates (see text for details), 
(middle) Decomposition of the same transform based on the 
gate U (0), designated with two joined circles. Note that U{cf)) 
can be achieved with simple free evolution, driven by the 
Hamiltonian (1). T(0) is a one-qubit phase gate. See Eqs. 
(7). In both cases not shown are swap gates at the end of the 
circuit, which exchange qubits 1 4 and 2 3. (bottom) 

Simultaneous conditional evolution is employed. 


where Rk{d) = Rk{d,0)- Sometimes we will drop the 
phase dependence in the gates, e.g. Uki{4>ki) will be writ¬ 
ten more concisely as Uki- 

We derive our implementation from the circuit shown 
in Fig. 1 (top), where we replace the CPHASE gates 
with free evolution U{(j)) and phase gates T{cj)). To 
this end, we note that a CPHASE gate Cki{(f>) act¬ 
ing on qubits k and I can be written as Cki{(p) = 
gi0|ii)fci(ii|fci ^ g*<^(|i)(i|)fc®(|i)(i|)i and, using the identity 

|I)(I| = (1 - cr^)/2, we get CkiiP) = 

After expanding the product we arrive at 

Cki{P)=ciTk{i)Ti{i)Uki{i), (8) 

where all gates commute. Hence, up to an unimportant 
global phase, which will be omitted in what follows, a 
CPHASE gate can be achieved with two phase gates and 
the free evolution gate. 
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Next we proceed with reconstructing the gate sequence 
from Fig. 1 (top), where our aim is to obtain a circuit 
involving only free evolution gates and a minimum num¬ 
ber of single-qubit gates. We note that the gate series, 
encircled with a dash line, can be expressed as 


Cl2(<('2)C'l3(03)C'l4(</>4) — 


exp 


Ti (^) T2 (^) T 3 (%) r4 (^) , 


(9) 


where ipi = 4>2 + (t>3 + 4’i or 4'i = St 2 ^or any N. If 
we substitute the phases from Eq. (5), we get 


<^i = f(l-2 -'^+i), (10a) 

(j)i = -K for I > 1. (10b) 


We concatenate such sequences of CPHASE gates with 
Hadamard gates, according to Fig. 1 (top), and then we 
rearrange the so-obtained circuit such that all phase gates 
to the left (right) of the Hadamard gates are combined 
in a single phase gate placed at the beginning (end) of 
the circuit. Thus we reconstruct the entire circuit in the 
form shown in Fig. 1 (middle). Therefore, the Fourier 
transform over iV = 4 qubits can be achieved with the 
following gate sequence (in operator form, to be read 
from right to left) 


Tf Ha e 4V2 


-( 1 ,(2)^(3)+ 1 ,( 2 ) ^(4)) 


li?3 e 
1 


j2:/'l,(l),(2),i,(l)^(3), 1 (1) (4)\ 

xH2e^ ‘ ‘ JHiTi, ( 11 ) 


or, in the general-TV case, 

Tf Hn Pn-1 Hn-1 ...P 2 H 2 P 1 Tj, (12) 


where 


,- 2 L 

= e* 4 ^ 1 =.. 


(13) 


We will refer to Eq. (12) as the consecutive sequence. 
Here is the Hadamard gate, applied to qubit k, and Tj 
and Tf are phase gates, applied to all N qubits at the be¬ 
ginning and at the end of the circuit, Tj = 

Tf = nf=i7’fc(</>F.fc), with = f(l - 2-'=+!) and 
'('F.fc = f (1 — In a more compact form we have 


• N 2L/1 — ^ \ 

= g*Efc=i 4(1-2 

E N ^ /i o —fc + l\ _(A:) 


(14a) 

(14b) 


Note that a phase gate represents a mere shift of the 
phase of the driving field, rather than a physical modifi¬ 
cation of the qubit. As such, it is not associated to any 
infidelity. 

In the rest of this section we show that by a proper tai¬ 
loring of the coupling matrix J it is possible to execute 
simultaneously the whole sequence Pk (cf. Eq. (13)) of 


conditional gates enclosed between two Hadamard gates 
(cf. Fig. 1). The tailored evolution operator U we de¬ 
note as U{J) in order to highlight its dependence on the 
coupling matrix J, which will different for different 
In units where t = 1, U{J) is defined completely by the 
matrix J. Now in place of the sequence (12) we use the 
following 

Tf Hn Hn-i ...t/(j(')) i?2 C/(J«) Ti, 

(15) 

which we call parallel sequence. 

Below we will determine the elements For con¬ 

ciseness we consider the implementation with four ions 
[N = 4), shown in Fig. 1 (middle). must satisfy the 
phase relations listed in Eq. (5): 

—_ '’L _ ('15') 

-'fei - 2i-k+i^ 

where fc = 1,..., A^ — 1 and k + 1 < I < N with iV = 4. 
Note that no net phase shift must take place between 
ions 2, 3 and 4 in the application of 17(J*-^^) and 3 and 4 
in the application of This leads to 

^ 3 ’ = Aa = 0, (17a) 

jW + = 0. (17b) 


In Eq. (17b) we account for the fact that 
and C/ 34 (J^^^) commute with the Hadamard gate on ion 
2. Therefore they can be combined into a single gate 
which must produce no net effect. 

It is important to note that the righthand side of Eqs. 
(16) and (17) is defined up to the addition of 27rz, where z 
is an arbitrary integer, z G Z. This offers larger flexibility 
for the values of J. 

Now we assume that can be factorized, = 
TTaa^^ a4\ where we introduce the physical parameter 
a in order to obtain dimensionless coefficients . Then 
a solution to Eqs. (16) and (17) for any N is 


4"’ = 0 

^(n) _ __ 

” ^{2qn - 9 „- i ) 2^-"+3 

- . 2“^+"“^-I-20*-"^ 

(™) _ f_ 4- ZCj. 


for k < n, 
for k = n, 

for k > n, 


(18a) 

(18b) 

(18c) 


(n) 

where qk and c), are arbitrary integers, 2qk ^ qk-i and 
qo = 0. A simple example is when qk = 1 (k > 0) and 



III. PHYSICAL REALIZATION WITH ION 
TRAPS 

Below we illustrate our method with homogeneous 
magnetic-gradient ion traps [16], where the implemen¬ 
tation of the [/-gate is appealingly simple: it occurs nat¬ 
urally in the time evolution of the energy eigenstates of 
the Hamiltonian (1). 
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A. Hamiltonian of the system and implementation 
of the gate set 

Consider a system of N ions, confined strongly in two 
orthogonal directions (x and y) and by a weaker potential 
along a third direction (z), known as the trap axis. The 
ions are laser cooled so that the ions form a linear chain. 
Two internal levels of each ion, which we denote as |t) 
and 14,), form a qubit. These can be a pair of hyperfine 
ground states, serving as effective spin-1/2 states. In 
addition, a magnetic field B = B{z)ez is applied, such 
that the ions experience a gradient along the trap axis. 


for example, when the spins are prepared in a particular 
state as part of the initialization of the system. Because 
the interactions that we will consider below do not couple 
the vibrational motion to the spin states, the latter are 
evolved independently of the motional states and thus 
p{t) preserves its product form over time. Therefore with 
U{t) we will refer only to the spin-driving part of the 
propagator, as we will be only interested in the dynamics 
of the logical basis states. For this reason we will omit 
the second term from the Hamiltonian (19). Moreover, 
because U {t) is independent of the vibrational state and 
it does not matter whether the vibrational state is cooled 
to the ground state or not. 


1. Hamiltonian of the system 


2. Implementation of the gate set 


The Hamiltonian of the system of N ions, up to second 
order in the ions’ vibrational motion, is given by 

, N N k ^ 

n—1 n—1 n<.m—l 

(19) 

The first term is the sum of the internal energies of the 
two-level ions located at equilibrium positions and rep¬ 
resented by Pauli’s matrices Due to the magnetic 

field gradient, the ions’ transition resonant frequencies 
are position-dependent, Wn = gp,BB{zn)/fi, where g and 
Pb are respectively the electron’s g-factor and the Bohr 
magneton. The second term describes N axial collective 
vibrational modes with frequencies and a„ are 

respectively the creation and annihilation operators of 
mode n. The last term realizes a long-range pairwise in¬ 
teraction between the spins, known as spin-spin coupling 
[16], whereby spin n is coupled to spin m with strength 


In our implementation we will use two bases, shown in 
Fig. 2. The coupled basis is formed by the magnetic-field 
sensitive states | -I- 1) and |0'). Ions residing in this basis 
are coupled by the magnetic-field gradient and thus can 
participate in conditional gates. The uncoupled basis is 
formed by the magnetic-field insensitive states |0) and 
I O'). We will uncouple ions when we want to suppress 
interaction with the others. Below we describe a possible 
scheme to implement coupling and uncoupling. 

To illustrate the scheme we consider the state vector of 
the three states \ip) = (coq c+i, co)"*". The general initial 
(coupled state) is IV'i) = (cos0, sin0,0)"^. We imple¬ 
ment a TT-pulse on the transition |0)„ O |0')„, 


IIo'.o = 


0 0 1 
0 1 0 
1 0 0 


( 21 ) 


and a 7r-pulse on the transition | -|- 1)„ o |0')„, 


Jn 


jgi^Bf dB{z) 
2h dz 


dB{z) 


dz 


{A-^) 


( 20 ) 


Here Amn is the Hessian of the total potential energy, 
which is the matrix of second-order partial derivatives of 
the external plus Coulomb potential taken at the equilib¬ 
rium positions Zn- Thus the spin-spin coupling is defined 
by the values of the magnetic field gradient and the shape 
of the potential surface around the equilibrium positions 
of the ions. In the following we will discover that only the 
third term is essential for the implementation of QFT. 

For the complete description of the system, we in¬ 
troduce the density operator pit), which in the ab¬ 
sence of decoherence, as assumed below, obeys the Li- 
ouville equation dp{t)/dt = —i[Hsit),p{t)], with Hs{t) 
being the Hamiltonian. The formal solution is p{t) = 
U{t)piO)U^{t), where 17(t) is the propagator (the evo¬ 
lution matrix). For our time-independent Hamiltonian 
(19) we have 1/(1) = exTp{—iHst/h). 

We assume that initially the system resides in a pure 
state, which is a product of a spin part and a motional 
part, p(0) = Pintemai(O) 0 ^motional (0). This happens. 


IIo'.+i = 


0 1 0 
1 0 0 
0 0 1 


Thus we obtain 


( 22 ) 


|'^)uncoupled — iio'.o n 0',-|-l nO',o|V’)coupled- (23) 

The second Hoqo-pulse automatically takes care of a 
potential problem that the Ho'^o-transition is Zeeman 
insensitive and thus swaps populations for all ions. It 
can be omitted, for speed-up and in the sense of ’fewer 
gates=lower errors’, but then some book-keeping is re¬ 
quired, where all the populations are. If applied, a 
useful improvement is then to do the second rotation 
|0)„ O |0')n with a phase of tt, so pulse length errors 
would be compensated as the Rabi frequencies for all 
ions on that transition might not be exactly identical. 

The U-gate. In a frame rotating with the spins’ 
transition frequencies our system is described by 
the Hamiltonian H'^ = F'^HsF — iHF^F, where F = 

exp{—it/2J2n=i^nO'^^'^) [19]. The transformed Hamil¬ 
tonian Hg contains only the spin-spin coupling term and 











5 


+ 1 ) 



FIG. 2: Two bases of quantum states are used. The cou¬ 
pled basis comprises |0') and the magnetic field sensitive state 
I -I- 1). The uncoupled basis encodes the qubits and is formed 
by magnetic field insensitive states |0) and |0'), which belong 
to different hyperfine manifolds. An ion is brought to the un¬ 
coupled basis to suppress its interaction with the other ions. 

thus it coincides with the Hamiltonian (1). We remind 
that so far no external interaction is applied to the ion 
chain other than the trapping electric potentials and the 
magnetic gradient field. Hence the entangling gate Uki {(j)) 
(cf. Eq. (7c)) occurs from free evolution of ions k and I 
for time tki = 2(j)/ Jki with all other ions being uncoupled. 

The R-gate. The one-qubit rotation Rk{0,4>) is ob¬ 
tained using a microwave (mw) pulse resonant with the 
spin-flip transition frequency of ion k. The mw pulse 
delivers an oscillating magnetic field in the x-y plane 

B' = B' [cos {ut + (j)) ex + sin (wt -|- </>) Cy]. (24) 

It couples to the spin thereby driving the rotation 

Rk {0,4>) = exp[— \i6 cos 4> + sin ], (25) 

where 9 = Qt is the pulse area and H = gyL^B' /h is 
the coupling between the spin and the driving field. The 
angle (j) is set by the phase of this field. 

The H-gate. The Hadamard gate, up to an unimpor¬ 
tant phase, is obtained with two pulses, 

= i?(7r/2,-7r/2)i7(7r,0), (26) 

where the second pulse is phased with —7r/2 relative to 
the first. 

Armed with a gate set, in the following we show how 
to realize QFT. 

B. Realization of the Fourier transform 

1. The consecutive sequence 

For brevity, we illustrate the implementation with four 
ions. The gate sequence is shown in Eq. (11) and is 
depicted in Fig. 1 (middle). After the application of 


the T-gates we uncouple all ions and follow the steps: i) 
Couple ion 1 and then apply the Hadamard gate H to it. 
Next continue with the implementation of the series Pi. 
ii) Couple ion 2. Now ions 1 and 2 can interact. Wait for 
time T = ( 7 r/ 4 )/Ji 2 thereby applying the gate C/i 2 ( 7 r/ 4 ) 
and then uncouple ion 2. Continue likewise with the 
gates 17 i 3 and t/ 14 . iii) After completing Pi uncouple ion 
1. Thus it is no longer manipulated until the application 
of the final phase gate (cf. Fig. 1 (middle)). Steps i)- 
ii) complete the encircled part of the circuit. Likewise 
proceed with the rest. In the general-case one follows 
Eq. (12). 


2. The parallel sequence 

We now discuss the sequence (15) depicted in Fig. 1 
(bottom), where we focus on the realization of the tai¬ 
lored operator U{J). This operator can be achieved using 
oscillating laser or magnetic fields, for example. Below 
we employ magnetic fields, while we note that with laser 
fields the implementation is essentially the same. 

Consider N trapped ions, as described above, which 
now interact with a bichromatic magnetic field with fre¬ 
quencies bJr = Wo,n — Wc.m. + d (“red”) and Wf, = Wo.n + 
Wc.m. — 5 (“blue”) (wc.in. ^ |<^|) timed near the centre-of- 
mass (c.m.) mode Wc.m. with detuning 5 (wc.m. |i5|)- 
Here wo,n is the atomic transition frequency for ion n. 
The interaction Hamiltonian in the Lamb-Dicke limit and 
in the rotating-wave approximation is given by [20-22] 
(^=1) 

N 

Hi = J2 5fca((^+)(ate*‘5*-*^'^ + ). (27) 

k=l 

Here gk is the individual (time-independent) spin-phonon 
coupling for ion fc, a) and a correspond to the c.m. vibra¬ 
tional mode and (T{Tk ) — o'fe with be¬ 

ing the spin raising and lowering operators for ion k. The 
spin and motional phases are defined as = ^{<fk + Tk) 
and = ^{fk ~ Tk)’ '"^here the phases and ifl. cor¬ 
respond to the blue and the red component of the oscil¬ 
lating magnetic field. Hereafter we assume that the blue 
and the red B-field are counterpropagating, which implies 
that = 0 and = pk, such that cr{p'^) = a^. 

The propagator is obtained using the Magnus expan¬ 
sion [23]: 

U = D{a) exp ’ (^8) 

where D{a) = and a = (1 — 9k^k- 

A more detailed study of interaction with bichromatic 
fields can be found in Ref. [20]. At time r = 2tt/5 the 
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displacement D{a) vanishes and we obtain 


U = exp 




exp 


N 

■ 9k9p<Jk(Tp 

k<.p—l 


(29) 

This propagator has the form as shown in Eq. ( 6 ) with 
Jki oc gkQi- The global phase factor can be easily com¬ 
pensated for at the end of the circuit. 


For homogeneous traps having nearly harmonic elec¬ 
trostatic potential and magnetic field gradient b = 20 
T/m we have Ti « 3.241 ms, T 2 ~ 0.558 ms and T 3 « 
4.478 ms. Due to the simultaneous interactions, the total 
duration of the algorithm is now Ttot = T 1 +T 2 +T 3 fv 8.3 
ms. For comparison, the duration using sequence (12) is 
tt/( 4 J 12 ) -I- 7r/(8 J 13 ) -I- tt/( 4 J 23 ) ~ 11-2 ms. For the pulse 
areas we have Ai « 0.6547r and A 2 ~ O.TTItt. 


C. Example with N = 3 qubits 


IV. DISCUSSION AND CONCLUSIONS 


Here we specify our implementation (12) for a system 
of three qubits. Following our proposal (cf. Eq. (12)), 
the three-qubit Fourier transform is obtained with 


The total implementation time for the consecutive se¬ 
quence from Eq. (12) is estimated by summing up the 
time it takes to produce each {7-gate, 


TFi73e”v8‘""'^" 


7726 


/I z z , 1 z ^ N 

ITT g CTj 0-2 -I- ig 0-1 0-3 j ^ 


(30) 


We can rearrange this sequence so that we can employ si¬ 
multaneous interactions in the first two conditional gates 
and thereby gain some speed-up. The rearranged se¬ 
quence is 

R2{A2, f )i?3(f, -f )t723(r3)i?i(7r, ff )i72(Ai, ^) (31) 
i?3(7r, - f|){7(T2)i?3(^)C/(Ti)i?i (7r)i?2(^)i?3(7r)7{-i, 

The derivation of this circuit is shown in Ref. [24]. Here 
we have used Eq. (26) and we have 

ri = I (i + ifc) . (32a) 

=f(i-dr;)^ (32i>) 

The duration T 3 and the pulse areas Ai and A 2 are ob¬ 
tained from the following set of equations 


1 fd-Tsla+A p.iJ23T3/2 

\pi 

sin ^ sin 



cos ^ cos 2 ^ = 0 , 

(33a) 


sin cos 2 ^ 6* ■^2373 


- 

cos sin 2 ^ = 0 . 

(33b) 


^-1 ^ n-(l-k+l) 

k—1 l—k-\-l 


Jkl 


(34) 


We obtain T{N) Ri cW^, where, for example, c is on 
the order of Milliseconds for experiments with 
gradients of ti Ri 20 T/m and an axial trap frequency 
of t'o Ri 27r • 200 kHz. The implementation needs 2{N- 
1) phase gates, N Hadamard gates and 0{N'^) 7 r-pulses 
for coupling and uncoupling ions, amounting to O(iV^) 
pulses overall. 

The time duration of the parallel sequence (15) where 
now the conditional interactions run simultaneously is 
given by 


AT-l 

= I E 75^- (35) 

n—1 ^n,n+l 

Using Eq. (18) we obtain T{N) = ^ v. Note that 

a quadratic speed-up is achieved. 
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